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There exist several theoretical motivations for primordial correlation functions (such as the power 
spectrum) to contain oscillations as a logarithmic function of comoving momentum k. While these 
features are commonly searched for in fc-space, an alternative is to use angular space; that is, 
search for correlations between the directional vectors of observation. We develop tools to efficiently 
compute the angular correlations based on a stationary phase approximation and examine several 
example oscillations in the primordial power spectrum, bispectrum, and trispectrum. We find that 
logarithmically-periodic oscillations are essentially featureless and therefore difficult to detect using 
the standard correlator, though others might be feasible. 

PACS numbers: 98.80.Cq, 98.80.-k, 98.70.Vc 



INTRODUCTION 

The leading paradigm in modern cosmology is infla- 
tionary theory, in which quantum field vacuum energy 
causes the Universe to undergo a very rapid expansion 
in a very short amount of time ^0]. One of the con- 
sequences of inflation is that the quantum field fluctu- 
ations source structure formations, which then become 
fluctuations in the CMB temperature and polarization 
as well as large scale structure. Thus, precision cosmo- 
logical measurements will yield tight constraints on the 
quantum field interactions, and hence the microscopic 
details of the inflationary theory. Such measurements 
are now being made with increasing precision by COBE 
[5] and WMAP 0(7], and there will soon be data from 
the Planck [H] and Euclid [9] satellites, and suborbital or 
ground-based polarization-dedicated experiments |10) . 

The simplest inflation model, and the one considered 
here, contains a single field <p which can be decomposed 
into the background and fluctuations as 

0(x,r) = 0o (t) + if(x,T). 

We use the conformal time — oo < r < CP such 
that t — y — represents future infinity, and then Fourier 
transform the fluctuations into comoving momentum k- 
space modes (fkir). The interactions of the underlying 
inflationary quantum field theory are then encoded in 
the correlation functions of (p^. It is these microscopic 
interactions which determine the quantum field theory 
model, and hence the precise mechanism responsible for 
cosmological inflation. While a variety of inflationary 
models have been proposed, differentiation between them 
requires precision measurements of possible features in 
the correlation functions. Such a measurement is then 
tantamount to explaining how the Universe came into 
being. 



While some of the correlations may contain features 
slowly varying as a function of k, others could be sharper. 
One such feature would be oscillations. This modulation 
could occur through high-energy physics encoded in the 
choice of vacuum [11H35] , sharp turns in multifield infla- 
tion models [37H12], or quasi-periodic backgrounds [151 - 
135] . While most attention has focused on oscillations in 
the power spectrum, there has also been increasing inter- 
est in the bispectrum [IrjHlfj] and trispectrum [Mj . 

This developing understanding of the relationship be- 
tween fundamental physics and primordial oscillations 
motivates an efficient method to detect such features. 
Such searches typically have analyzed the correlation 
function as a function of /,• space 50 56 . This com- 
parison is usually done by transforming the primordial 
perturbation into the CMB temperature fluctuation at 
degree t. In the approximation of instantaneous CMB 
decoupling at time r\ the transformation is given by 

where T ra( j(fc) is the radiative transfer function. Here 
we have simplified the gauge-fixing and freeze-out issues, 
which are irrelevant to the discussion at hand. Each 
^-mode receives contributions from angular separation 
80 « ir/£. Correlations of ATg/T can then be evaluated 
using correlations of ip^, such as the power spectrum 




While this is indeed the natural way to constrain slight 
scale-dependence, it is not so natural for detecting oscil- 
lations. 



A more promising approach for detecting oscillations 
would be to Fourier transform to angular, or direction, 
space [57] via 

This has the advantage that the observation is completely 
localized to a single direction, and so one no longer needs 
to sample many region of the sky to gain information for 
a single £-mode. This may also produce a stronger sig- 
nal and allow a more efficient computational processing. 
Such an approach has not yet been widely explored, and 
is well worth investigating. 

In this article we calculate the angular correlation func- 
tions for models containing oscillations. We focus upon 
oscillations which are periodic in the logarithm of k since 
these are physically motivated by new physics arising at 
energy scale M. We use a stationary phase approxima- 
tion which takes advantage of the oscillatory features. 
From this we derive the dominant configurations of the 
direction-vectors and explicitly evaluate the correlation 
function at these locations. The result is correct to lead- 
ing order in H/M, where H is the Hubble parameter de- 
termining the energy scale of inflation. To maximize the 
possible benefit of our new approach we assume instan- 
taneous decoupling; were we to relax this assumption, 
oscillations more rapid than the timescale of decoupling 
would become smeared out. In §1 we evaluate an ex- 
ample in the primordial power spectrum. In §2 we apply 
this technique to an example in the bispectrum, and in §3 
we evaluate the trispectrum. In §4 we discuss the results 
and offer observational prospects. 



2 





\ A AAA A A A A A A A A A A A/ 




\ 1 \ 1 W l ) \ 1 11 If i 

Vl/VVVVl/V 


V V V 

2/rA 


\ 1 II II 11 

\l V V V V 




k 



FIG. 1: Primordial power spectrum containing oscillations 
periodic in lnfe. For reference, the scale-invariant "flat" power 
spectrum is given by the thin red line. 
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FIG. 2: Simplified CMB temperature power spectrum in the 
^-basis. The flat power spectrum is given by the thin red line. 



POWER SPECTRUM 



fluctuation power spectrum, 



Let us begin with the power spectrum, or two-point 
correlation function: 



<*k(0-)M0-)) = ^^W(2^) 3 ( 5 3 (k + k'). 

Suppose the power spectrum is mostly scale-invariant 
but contains oscillations set by the scale H/M, which 
is strongly theoretically motivated [TTH3S]: 



P v (k)=P 



. M k 
1 + p cos — In — 
H k 



(2) 



where fco is some reference wavenumber and /3 is a 
small parameter controlling the oscillation magnitude, as 
shown in Figure [I] It was shown in |36j that there is a 
theoretical bound on the oscillation frequency. 

A quantity of great interest is the CMB temperature 



C, 



(2tt) 3 (2tt) 3 



(3) 

{cpktpv}T lad (k)ji(r]k)T r:id (k')ji(rik'). 



The transfer function T ra( j(fc) is in general a complicated 
function of scale and cosmological parameters, but at 
both large and small angular scales it can be approxi- 
mated as T ra d(fc) e~ k l k ° for some effective diffusion 
scale kjj- For simplicity of an analytic answer we will 
use this function, and the acoustic peaks appearing at 
intermediate scales would be straightforward to include 
numerically. Evaluating this with the primordial oscil- 
lation ([2]) gives the power spectrum shown in Figure [2j 
demonstrating that oscillations now appear as (damped) 
oscillations in the Cg. 

An alternative, though much less common, represen- 
tation is to evaluate the temperature fluctuation power 
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spectrum in position space, 
AT(ni) AT(n 2 ). 



(A 77 7)2 



T T 
d 3 k d 3 k' 

(2tt) 3 (2tt) 3 



(^k<Pk')T 1 rad (A:)e 



1 Trad(fc')e 



-tJ7k'-n2 



We now define the rescaled conformal time 7 = Hrj/M, 
and use the unit vector difference 

N = n! — n 2 

so < N < 2. Due to isotropy, the power spectrum can 
only depend on N. Writing the cosine as exponentials, 
the correction term can now be cast into the form 



/ AT V 
where 



d 3 k 

(2tt) 3 



G(k) 



-if F_(k) 



G(k) 
F±(k) 



fc 3 

fc 

± In — + 7k • N. 

fc 



The majority of the contribution will come from near the 
extremum at 



fc*± 



This is easily solved to yield 



k*± = =F 



N 

^N2- 



At this moment of stationary phase, 

G(k»±) = ir^Poi-fNfe- 2 ^^ 2 , 
F±(K±) = T[ln( 7 ^ ) + 1] 

and the fluctuations around this point are given by 

d 2 F ± 



Mf 



dk a dk b 



k, 

j2 cab 



±7 2 (N 2 5 ab - 2N a N b ) 



After performing the Gaussian integral the answer is 
then conveniently written as 



/AT 
(~T~ 



d 3 k 



(2tt) 3 



(iH) 3 



{2n) 3 M 3 \M 




FIG. 3: Simplified CMB temperature power spectrum in the 
direction-basis, as a function of directional separation N. The 
flat power spectrum is given by the thin red line. 



To explicitly evaluate this we use rotational invariance to 
choose coordinates such that 

N = Ni 

which makes 



\M±\ = ( 7 A0 e 



The final answer for the oscillation term is then 



m - 

\ / osc 



= 2tt 2 /3P 



H 

2ttM 



3/2 



M 

cos I — [In (7.^0) 



(4) 



+ 



This is shown in Figure [3j Contrary to the expectation 
that there would be a single large peak, there are again 
multiple oscillations. This is due to the primordial power 
spectrum being logarithmically periodic, which dilutes 
the signal over all scales. We will find this behavior con- 
tinue to higher-point correlators. 



BISPECTRUM 

We now consider a more elaborate calculation, the 
position-space primordial bispectrum: 

AT(ni) AT(n 2 ) AT(n 3 ) 

\ rj~\ rjn rjn I 

r 3 d 3 ki 

= J II J^j3 T ™ d ( k ^ eX P ( _ "? k * ' A i) (VkiWkaVO 

where the correlation is of the form 

(PkxPkaVO = S v (k 1 ,k 2 ,k 3 )(2^) 3 ( 5 3 [J2*i J • 
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Since ki + k 2 + k3 = the fc-space correlations are cat- 
egorized by the triangle formed by the k^. Oscillations 
(or other features) could then depend arbitrarily on both 
the size and shape of the triangle. Let us consider the 
motivation for these two types of bispectrum oscillations. 

The primordial fluctuations have been found to be 
very nearly scale- invariant. Though oscillations in the 
power spectrum (by definition) break this invariance, it 
is possible to maintain perfect scale-independence and 
yet have oscillations in the bispectrum. This could be 
accomplished via a phase which depends only upon the 
relative angles formed by the k^ and not by their overall 
magnitude, being invariant under fcj — > Xki. Despite this 
strong theoretical motivation, it appears difficult to for- 
mulate any non-trivial scale-independent bispectrum os- 
cillations which allow a stationary phase approximation. 
Thus here we will limit ourselves to scale-dependent oscil- 
lations. These would contribute to other scale-dependent 
features in the bispectrum, and which are possibly mea- 
surable in the near future 1551. 



The extremum will satisfy 



dkf k 2 



K *±,i 

1 + ^*±,2 + ^*±,3 



+ 771° + W a 



For now we will treat w as a constant and solve for \c*±a 
in terms of it. This yields 



k*± a — ^ 



7nj + w 



Ei h^i + w| 

The ki-fluctuation matrix near this point is 



(8) 



M 



k,ab _ u r ± 

±>y ~ dkfdk b 4 

1 3 




| 7 n fe + w| 2 6ijS ab T 2 (-yh? + w a ) (7^ + w b ) . 



Scale-Dependent Oscillations 



Substituting the solution ^ back into the exponent 
F± in yields 



In slow-roll inflation the bispectrum induced by grav- 
itational interactions is negligibly small [5§| so we can 
focus exclusively on the term generated by field interac- 
tions of a particular model. Consider the example 



B- 



Bn 



(hk 2 k 3 y 



k\ 



k 2 



k 3 



3/cq 



c.c. 



(5) 



This is similar to the "Resonant Non-Gaussianity" model 
studied in |48j . As in the previous power spectrum ex- 
ample, the (kik 2 k 3 )~ 2 compensates for the phase-space 
factor; the bispectrum ^ then depends only on the sum 
E» k 2 , and so is well-behaved for any non- vanishing tri- 
angle. 

In order to perform the integrals over kj we first employ 
the Fourier representation of the momentum-conserving 
delta- function, where we include a factor of M/H for 
future convenience: 

{2,f8* (jy^ = (f)7 d 3 wexp _,| 

(6) 

Performing the same rescaling as before, the functions 
are now 



By 



{kik 2 k 3 y 



(7) 



F ± (k i)W ) = ±3 In 1 3fc 2 2 3 



i=l 



E k 0- 



F±(w) = t 



■In 3fc 2 ]T| 7 n 4 +w| 2 



We then treat w as another parameter to vary, whose 
extremum yields the momentum conservation equation 



_ 9F± _ Ej=i 7"? + «tf 



dw a 



E,- 



2 ' 



(9) 



Note that we could have obtained this by taking 
dF±/dw a in ([T]) originally, but obtaining it in this order 
makes evaluation of the fluctuation-matrix determinant 
much easier. The solution to (|9| is 

7 \ - - 



The fluctuations around this point are 
d 2 F ± 



m: 



dw a dw b 



±7" 



tab 



The fact that M™ ~ (A4 k ) fits our intuition that w 
is a "negative degree of freedom," constraining the kj to 
be conserved. 
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The temperature bispectrum is then 




FIG. 4: The directional vector configuration (top) resulting 
A Family of Solutions m a fe-space equilateral triangle (bottom). 



While analyzing the general bispectrum solution ( 10 ) 
is difficult, it is simple to consider a one-parameter family 
in which the symmetric directional vectors all make an 
angle 9 with some z-axis as shown in Figure [4j 

27TZ 27TZ 

Aj = sin 9 cos x + sin 9 sin y + cos 9z. 

3 3 

This makes the solution 



w* 


= —7 cosy, 




(cos ^px + sin ^y 
37 sin 9 


J&i + | 


= 7 sin 9, 


*±(w.) 


= =F [ln(7fc sin6 l ) + 1] . 




= ±(3 1 2 sm 2 9)~ 1 6 ab 1 


\M^ ab \ 


= (37 2 sin 2 6>)- 3 , 




= ± 7 2 sin 2 9 [3%<5° 6 



2iri „ . 2ni „ 
2 j cos — ^— x + sin ~^~y 



IM^fA = (3 7 2 sin 2 0) 9 . 



27TJ „ . 27TJ A 

cos x + sin y 

3 3 y 



The final answer is then very reminiscent of the power 
spectrum result Q but with N — > sinf?, 



/AT 



' \2irM J 
M 

■ vi i> 1 — [ln(7fc sin6') + 1] 
H 



TRISPECTRUM 

The trispectrum is not as well-studied as the bispec- 
trum but should soon also be an important tool in ana- 
lyzing models of inflation [60ti63j . The four-point corre- 
lation in position-space is 



,Ar(ni) AT(n 2 ) AT(n 3 ) AT(n 4 ), 

\ rj~i rj~i rj~i rj~i 1 



(11) 



= j n (2^ Trad ^ fe ^ exp (~^ ki ■ (^k! vk 2 ^k 3 ^k 4 ) 

where 

(yki^kaPksPkJ = 7 1 v (k 1 ,k2,k 3 ,k 4 )(27r) 3 (5 3 ^]C ki ^ ' 

The only theoretically motivated trispectrum oscillations 
were derived in 34 , in which a light field in an excited 
state interacts with a heavy field, as shown in Figure [5] 
(here we are simplifying a bit to avoid subtleties of the 
"in-in" formalism). The correlation is given by 

T = 9 x (12) 

v 2 10 7r 2 (fci/c 2 fc 3 fc4) 2 M 1 ; 

[(h + k 2 f |k x + k 2 | 2 ] ~ 1/4 [(k 3 + k,f - |k 3 + k 4 | 2 ] " 1/4 x 



So the bispectrum will also be suppressed at small 



h + fc 2 + y/(fc! + fc 2 ) 2 - |k! +gg 
fe 3 + fc 4 + V / (^3 + fc4) 2 - |k 3 + k 4 | 2 / 

+permutations. 



-iM/H 



C.C. 
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The First Half: Solving for ki, k 2 

Taking the derivatives of Fl with respect to ki, k2 re- 
sults in the following extrema conditions: 



fci 



/|'4 



ki 



fe 2 (k 1 -k 2 ) 



y / 2fe 1 fe 2 (i-k 1 -k 2 ) 



FIG. 5: Tree-level high-energy contribution to the 4-point cor- 
relation function. Single lines denote light fields, the double 
line denotes a heavy field. The total correlation includes the 
two additional permutations of ki, k 2 , k3, k4. 



This is scale invariant under fcj — > Xki. Let us focus on 
the permutation written out, since the others are simple 
variations of these. 

In order to solve the stationary phase constraints, note 
that (aside from momentum-conservation) in the trispec- 
trum correlation ( 12 ) it is only ki and k2 which are cou- 
pled only to each other, and likewise for k3 and k 4 : 



2fcifc 2 (l-ki-ki 



k 2 + 



fci(k 2 -k!) 



l-ki-ka 



fcl 



2fc 1 fc 2 (i-k 1 -k 2 ) 



7&i + w, (14) 



7n 2 + w. (15) 



In analyzing these equations, it will be helpful to make 
a change of variables: 



K = ki - k 2 



N = ni — n 2 , 



v = — w + ni + n 2 . 
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T v (2nf6 3 J> = 



(13) Rewriting (14) and (15) in these variables yields 



M 

n 

where 



J d 3 w [T^ 



k 2e -^i(k!+k 2 )-W 



r k 3 ,k4 e -^i(k 3 +k 4 )- 



ki 



fc 2 K 



Tkl ,k 2 _ glg 3/2 rrad(fcl)T rad (fc 2 ) 



2 5 7rVM(fcifc 2 ) 2 



[(fci + fc 2 ) 2 ~ |ki+k 2 

iM/H 



,21-1/4 



ki + k 2 + Ky/kjib 

k 2 - fc * K 

h+k 2 + Ky/hk^ 



(v + N) , 
(v-N) 



(16) 
(17) 



x fci + k 2 + yJ2kik 2 (l - ki • k 2 ; 



The T , ( k3,k4 component is identical except for the opposite 
oscillation phase. The factoring of ( 13 ) allows tremen- 
dous simplification in solving the stationary phase con- 
straints by expressing F as the sum of a "left" half con- 
taining only ki and k 2 , and a "right" half containing only 
k3 and k4, and finally also the prefactor G: 

F(ki,w) = F L (k 1 ,k 2 ,w) + J F 1 fl (k 3 ,k 4 ,w), 
G(ki) - G L (k 1 ,k 2 )G fl (k 3 ,k 4 ) 



where we have used K = y 2(1 — ki • k 2 ). Dotting (16) 
into (17) gives 

KN 

k 1 +k 2 + Ky/kik 2 2 



= 1 ( N 2 - v 2) 



(18) 



There are now two linear combinations of ( 16 ) , ( 17 1 which 
will prove useful. Subtracting, 



K 



Ky/k^k 2 



7 N. 



where 



.9^ 3/2 T rad (fc 1 )r rad (fc 2 



Squaring this gives 



Thus 



kik 2 



2 5 ny / M(k 1 k 2 ) 2 
( 



F L = In 



2fc!fc 2 (l - ki ■ k 2 ) 
(l - ^ • k 2 ) 



1/4 ■ 



K 



Now adding J16J) and (171, 



h + k 2 + J2hk 2 1 - ki • k 2 



V 



7 2 7V 2 . 



-N. 



ki + k 2 + Ky/kik 2 



7 v. 



(19) 



(20) 



(21) 



(22) 



+ 7 (ki • ni + k 2 • n 2 ) + (ki + k 2 ) ■ w 

and similarly for Gr and Fr with ki i2 —¥ ^ t 4, and a sign 
difference in the logarithmic term of Fr. 



Since (ki + k 2 ) ■ K = it is clear that v should be 
decomposed according to 

v = mM + nN 
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where M contains a U(l) symmetry parameterized by 
the angle <fr, 



M(0) 



, ni + n 2 

3S <h . 

ki +k 2 

V4 - K 2 ' 



sin 0ni x ri2 



(23) 



Using (18), the vector constraint (22) can then be de- 



composed into the component along M, 

, 2 2 2N 

IN — m — n I = 

2iOV V j 



and along N, 



7(fci - fc 2 ) 
2K 



(N 



2 2 2\ 

m — n 



(24) 



(25) 



Eq. ( 24 ) is then easily seen to be 
2N . „ 



NK 



V^-K 2 
2N 



sin ( 



cos 9. 



Vi-K 2 ' 



Note that because solutions to ( 18 ) exist only for v 2 = 
rn 2 +n 2 < N 2 , and (24) only for m > we limit ourselves 



to a certain interval of tp, 

^max - 



K 

arcsin — . 
2 



If we define 



y 



K 

V4-if 2 



this can be inverted as 

K{r) = 



2r 



The solution for v can then be written as 

v{r,ip,4>) = (N\/r 2 + 1 sin ip - Nr ) M(</>) (26) 



where the range is given by (note that r is not actually 
the radius, merely a convenient parameter) 



< 4> < 2ir, 
where 



<P<Pu 



< r < oo 




FIG. 6: Solutions to v. Each point in the sphere of radius iV 
corresponds to a unique solution of ki , k2 . The color indicates 
the time of interaction r. The vertical axis is the direction N 
and the transverse directions correspond to M. 



Thus the position-space trispectrum at decoupling ra- 
dius r\ really does project onto only those interactions 
taking place at |r| ~ r\. Figure [6] shows the space of 
solutions to v. 



It is then easy to obtain 



k x + k 2 
h - k 2 



K ( 2N 



1 



2Kn 



j(N 2 -v 2 )' 



Then linear combinations produce 



ki 
k2 



K 
2^ 
K 

2l 



2N 



N 2 



v 2 i 1 n) 



2N 



N 2 



n 
N 



1 

N 
1 

N 



(27) 



Using (21) and (23) we can also solve explicitly for the 



unit vectors, 



The solution to v thus fills a sphere of radius N. Each 
r represents an (American) football-shaped surface sus- 
pended between v = ±N corresponding to a contour of 
interaction time. This can be seen by 



M 



(K^k 



MjN -qN 



HK(r) K(r)' 



k! = — fi/4 - K(r) 2 M + K(r)N 



(28) 



k 2 = — ^4 - K (r) 2 M - K{r)N 



We can of course combine them to get ki , k2 but these 
are not so interesting by themselves. More important is 



the combined momentum, 

V4- K 2 ( 2KN 



ki + k 2 = -- 



2 7 \N 2 -v 2 ~~N 
K2n -N 



M 



7(7V2 _ ^; 2)- 
i\: 2 (jV 2 +^ ; 2 ) 2{KN) 2 n 



j(N 2 -v 2 ) 2 j(N 2 -v 2 ) 2 
4 



N 



7 [(iV 2 + w 2 ) 2 - 4n 2 N 2 } 

x [-(TV 2 + v 2 )v + 2NnN] (29) 

where we have used p4| , 

4(7V 2 - t- 2 ) 2 



X(r) z 



The fluctuations are 

iL,ab 



M 



{N 2 + v 2 ) 2 -An 2 N 2 ' 
d 2 F L 



6 7 



(30) 



which are given by 



pa | k 2 (kl-k%) 




lb | fe 2 (fc?-^) 


1 V 2 *-'iMl-ki-k 2 )_ 




1 V 2fe i fe 2(l-ki-k 2 )_ 



S ab -k'k° 



+ 



h + k 2 + V 2 fcife(l-ki-k 2 



V2fclfc2(l-k 1 -k 2 ) (2fe 1 fe 2 (l-k 1 -k 2 )) 



3/2 



fci +fe + V2feife 2 (l - ki -k 2 ; 



L.ab 
12 



fc 2 (fc?-fcg) 
-v/2*:ifc2(l-ki-k2) 



fc 2 



/c j ^ ^2 — ^ 1 ) 

V2fcife2(l-ki-k 2 ) 



fci + fc2 + \/2fciA: 2 (l-ki-k 2 



V2fci*=2(l-k 1 -k 2 ) (2fc 1 fe 2 (l-k 1 -k 2 )) 3 



fci + fc 2 + V 2fc i fc 2(l-ki-k 2 



where reminds us to evaluate the answer using the 



expressions (27), (28) in terms of v and N. Of course 



A^ 22 ab is analogous to VW 11 ' ah with ki k 2 . 

The Second Half: Solving for k3, k 4 

The identical logic carries through for k 3j4 , except 
there is an overall minus sign for the LHS of (14), (15) 



and so the solutions are simply the opposite of ki 2 : 
k 3 = - (V 4 - if(r 3 4) 2 M 3 4 + #(r 34 )N 34 ) , 
k 4 = \ (V4-^(r 34 ) 2 M 34 - K (r 34 )N 34 ) . 




FIG. 7: Solution to the auxiliary variable equation of motion 
which must satisfy both halves of the system. 



Here we have labeled N 34 = n 3 — n 4 , as opposed to the 
original N = Ni 2 , and used the appropriate coordinate 
r 34 . The magnitudes fc 3 ,/e 4 are analogous to &i,fc 2 , re- 
spectively. The fluctuation matrix A4 R is opposite to 
M L given in ([30]), 



M R,ab = _ M L,ab 

with k 3j4 replacing ki j2 . 

Solving the Complete System 

At this point we have two separate sets of solutions: 
one for ki 2 and one for k 3 4 . Substituting them back 
into the exponent function F and returning Vi 2 ,v 34 to 
the variable w gives 



F(w) 



In 



|w + 7A1I w + yn 2 | 



(31) 



|w + 7n 3 | |w + 7n 4 | 
Taking dF/dw a yields the stationary phase constraint 

7n" w° + 7n 2 u>" + 7^3 w% + jn% 







•7111 



- 7 n 2| 



7 n 3| 



w* + 7n 4 1 
(32) 

As in the bispectrum case, this constraint represents the 
conservation of momentum and could have been obtained 
by setting £\ k, = 0. 

We desire a solution to w which satisfies both the left 
and right halves, as shown in Figure[7] Note that because 
w must be a solution to both halves, it must correspond 
to a v 12 which must be both within the sphere of radius 
N12 centered at — (fti + n 2 ) as well as a v 34 within the 
sphere of radius 7V 34 centered at — (n 3 + n 4 ). 

To facilitate solution construction, a nice interpreta- 
tion of the system can be made in terms of (logarith- 
mic potential) electrostatics as follows. Consider four 
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charges, two positive and two negative, placed at loca- 
tions fii on the unit sphere. Then F represents the po- 
tential of a test charge at location — w/7 in this back- 



ground, and the condition (32) is the condition for the 



test charge to be in equilibrium, stable or unstable. 
The fluctuation matrix for the auxiliary parameter is 



M 



w,ab 



o 2 f 



dw a dw b 



i=l 



■ 7 n i 



[|w„ + yhtfS^ - 2« + 7 n?)K + 7^)] 



where q±.2 = +1 and 53,4 = — 1. Performing the Gaussian 
integral, the final answer is 



1 

permutations. 



y/\M L \ \M R \ \M^ 



c.c. 



Finally we must include the contribution of the flat power 
spectrum. This amounts to simply summing the permu- 
tations of the power spectrum-squared, 



■aussian 



(ATy /AT\ 2 



V T J Ni 

2 



(AT 



N 34 

(") ! 

V 1 /AT,. 



(AT 



N 2 , 



A Family of Solutions 



While a general solution to (32 ) appears to be difficult, 
there is one family of solutions which is obvious. In the 
language of the electrostatics analogy, place the two pos- 
itive charges on the unit sphere at rii j2 around the z-axis 
symmetrically at angle 9 so that they are separated by 
length N, and use their difference vector N 12 = n x — n 2 
to define x. Now place the negative charges at 63.4 at the 
same height, but rotated by some angle <fr with respect 
to x. This makes the directional unit vectors 



ni 
n 2 
n 3 
n 4 
N 12 
N 34 
Mia 



sin 9x + cos 9z. 

— sin 9x + cos 9z, 

sin 9(cos (fix + sin (fiy) + cos 9z, 

— sin 9 (cos (fix + sin (fiy) + cos 9z, 
2 sin 9x, 

2 sin #(cos (fix + sin (fiy ) , 
-M34 = z. 



(33) 



This is shown in Figure[5J Clearly a test charge located 
at any point on the z-axis will be in (unstable) equilib- 
rium, the force from the positive charges canceling that 




FIG. 8: A family of solutions to the trispectrum stationary 
phase constraints given by eq. ( 33 1 . 




k u kz 



FIG. 9: The fc-space configuration corresponding to the same 
family of solutions. 



from the negative ones. Thus we consider solutions of 
the form 

= w z z. 

The fact that v 2 < N 2 means the solution must be within 
the sphere center located at 

— (ni + n 2 ) = -(n 3 + n 4 ) = -2cos#z 
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and have the bounds 



sin 9 < 



cos 9 < sin ( 



Denoting 



a = sin 2 9 + (1///7 + cos 0) 2 



and A4 w as the fluctuation matrix without the z- 
component, we have the solution shown in Figure [9] 



|w + 711*1 

ni2 
TO12 



7 2 (7, 



n 3i = 0, 

— m 34 = 2 (w z /j + cos 6) 

w z 1 7+cos 



sin (9 



1 



K 



2 r 



2(u; 2 /7 + cos6') 
sin 2 9 — (w z /7 + cos ( 
K 2 ^ 



M' 
G(w*) 



= arccos 1 — 



16 sin 9 sin 



7 4 cr 4 



= 0, 



2g 2 H 3 ( 7 sin6») 9 e-(T fc ° sine ) 2 



(2tt) 2 MK 



1 



27 sin 



ki = 



k, = 




cos 0x + — sin 



X ^ X • ^ 
— cos0x - — sm0y 




The fluctuation matrix components are 



M L,ab 



*L,ab 



+ 7 tr 



(7 




k^ k^ 


fL,ba 


/ 7(7 


'21 


— 


\2sin9 


a 


(s ab 


— k^k^ 



x a x b 

KI K X X 



k^ x J I k<2 ~\~ x 

1 



K 



M 



L,ab 
22 



(s^s)' (*+*■) CM 



j 2 a 



' ™2 "-2 



if / X 



Since 



H,a6 



-M 



L,ab 



, , l . . , , the magnitudes of their deter- 

minants must be equal, jA^^I = 

Evaluation of the fluctuation matrix determinants is 
too involved analytically, and so we resort to numerical 
methods. Assembling everything together, and remem- 
bering that w z is unconstrained, the correlation will be 



/ AT\ / H 

\ / OSC \ ' 

r 7(sin0-cos0) ^ z 



(34) 



G 



7 (sin0+cose) 2vr ^/\M L \\M R \\M'™\ 

We should also include the contributions from permut- 
ing the ki in the interaction diagram shown in Figure [5] 
This is accomplished merely by permuting the directional 
vectors n^. Consider again the electrostatics analogy pre- 
sented earlier in which rii , n2 denote the location of pos- 
itive charges on the unit sphere whereas A3, A4 are the 
location of negative charges. For the family of solutions 
under consideration this yielded equilibrium points for a 
test charge located along the z-axis. If we now permute 
the charges so that ni, A3 are positive charges and A2, 
A4 are negative, simple electrostatic intuition shows that 
there will cease to be any equilibrium points (except in 
the special cases in which the n* overlap, which are sin- 
gular anyway). Thus there is no soution to the stationary 
phase equation, and the contribution vanishes to leading 
order. The same holds true for the final permutation of 
ni, A4 being positive charges and A2, 113 being negative. 



Thus we do not need to modify the answer ( 34 1 
The gaussian trispectrum result is given by 



AT 



1 J N= 



± / N=2sin6 
2 



JV=^/2(l-cos< 



AT 



2 



Af=^/2(l+cos <f>) sin 6 



In Figure 10 the trispectrum is shown as a function of 



9, keeping fixed 4> = n/2. In Figure 11 this is evaluated 
as a function of <p, keeping fixed 9 = tt/2. 



DISCUSSION 

We calculated the 2-, 3-, and 4-point angular correla- 
tion functions for models with logarithmic oscillations, 
and found that they are essentially featureless. Thus the 
direction-space representation does not appear to offer 
any advantages in detecting their presence compared to 
the conventional fc-space analysis. 

Consider instead the power spectrum arising from a 
boundary effective field theory [T2], meaning the oscilla- 
tions are periodic as a linear function of k, 



AP V = AkcosCk. 
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(A 77 T)* 




tt/4 nil 

FIG. 10: Simplified CMB trispectrum in direction-space, 
varying 6 while keeping fixed (j> = tt/2. The flat power spec- 
trum is given by the thin red line. 



(AT IT) 4 




tt/4 tt/2 



FIG. 11: Simplified CMB trispectrum in direction-space, 
varying <f> while keeping fixed 9 = tt/2. The flat power spec- 
trum is given by the thin red line. 



The position-space correlation would then peak when 
(^rj ~ B J d 2 k t S (C - rfa ■ n,) . 

The solution to this forces the to take the same relative 
configuration as the but rotated an overall angle away, 
reducing this to a one-parameter family of solutions to 
examine. 

Similar simplifications for the linear oscillations may 
be possible if modified transforms are used, such as a 
transformation from £-space back to angular scales with 
a logarithmic scale: 

t 

with Pi n i (cos 9) the Legendre polynomial. Oscillations 
which are periodic in the logarithm of k such as ([2| will 
then appear as peaks at 9 w N, 

P OS c(9)~l35 (Ov-^Y 

Such alternative search strategies merit further investi- 
gation. 
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Fourier transforming this yields peaked correlations lo- 
calized to definite angular separations, 



/AT 



AS (C - Nri) . 



These have the advantage of only requiring certain config- 
urations of angular separation, and thus are vastly more 
computationally efficient. This would also work well for 
the bouncing cosmology model studied in |43j , which ob- 
tains a complete power spectrum (not just a correction) 
of the form 



P„ = Ak r ' 



cos 



k 



<P0 



For higher-point correlations this is even more efficient. 
Suppose the bispectrum were of the form 



B, 



2 , 2,3 cos [C(ki + k 2 + k 3 )} 
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